Effect of the Zero-Mode Landau Level on Interlayer Magnetoresistance in Multilayer 

Massless Dirac Fermion Systems 
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We report on the experimental results of interlayer magnetoresistance in multilayer massless Dirac 
fermion system a-(BEDT-TTF)2l3 under hydrostatic pressure and its interpretation. We succeeded 
in detecting the zero-mode Landau level (n=0 Landau level) that is epected to appear at the contact 
points of Dirac cones in the magnetic field normal to the two-dimensional plane. The characteristic 
feature of zero-mode Landau carriers including the Zeeman effect is clearly seen in the interlayer 
magnet oresist ance . 

PACS numbers: 71.10.Pm, 72.15.Gd 



Novoselov et al. and Zhang et al. have experimen- 
tally demonstrated that graphene is a zero-gap system 
with massless Dirac particles [T, They discovered 
new steps of the quantum Hall effect attributable to a 
zero-gap energy structure with linear dispersion. How- 
ever, graphene is not a bulk crystal but a single layer 
of graphite. The bulk crystal zero-gap system was real- 
ized in the organic conductor a-(BEDT-TTF)2l3 under 
high pressure [1 i, i, S 0, Hi- The existence of the 
massless Dirac fermions state in bulk crystal has alr eady 
been discovered in graphite and Cdi_xHgxTe 

But graphite is not a zero-gap system but a semimetal. 
On the other hand, Cdi_xHgxTe is a three-dimensional 
(3D) semiconductor with narrow energy gap of more 
than 6 mcV. In this sense, a-(BEDT-TTF)2l3 is the first 
bulk 2D material with zero-gap energy bands. Thus, a- 
(BEDT-TTF)2l3 provides us with a testing ground for 
the transport of the multilayer massless Dirac fermions 
system. In this Letter, we demonstrate that the charac- 
teristic feature of transport in multilayer Dirac fermions 
systems is clearly seen in interlayer magnetoresistance. 

A crystal of a-(BEDT-TTF)2l3 consists of conductive 
layers of BEDT-TTF molecules and insulating layers of 
I3 anions [12]. Since the conductive layers are separated 
by the insulating layers, carriers in this system have a 
strong 2D nature. The ratio of the in-plane conductivity 
to the interlayer conductivity is about 10'^. Therefore, 
each conductive layer is expected to be almost indepen- 
dent. Under ambient pressure, this material behaves as a 
metal down to 135 K where it undergoes phase transition 
to a charge-ordered insulator. In the low temperature 
phase, a horizontal c harg e stripe pattern for -1-1 e and 
is formed [H, [13, [l^, [T6| . When the crystal is subjected 
to high hydrostatic pressure above 1.5 GPa, the metal- 
insulator transition is suppressed and the metallic region 
is preserved down to low temperatures [13, [H, [I^ ■ This 
change in electronic state is accompanied by the disap- 
pearance of the ch arg e ordering state, as shown by Ra- 
man experiments [lp|. Note that in the high pressure 
phase, the resistance is almost constant over the temper- 
ature range of 300 K to 2 K. The conductivity for each 
conductive layer is estimated to be nearly equal to e^//i. 



In contrast with the resistivity that is almost constant, 
carrier density and mobility exhibit strong temperature 
dependence. Carrier density (n) estimated by the Hall 
effect is proportional to the square of temperature (T) as 
n txT"^ over a wide temperature range. Such a peculiar 
character of this material remains an issue to be resolved. 

Recent band calculations by Kobayashi et al. sug- 
gested that a-(BEDT-TTF)2l3 under high pressure is in 
a zero-gap state 0, Q- The bottom of the conduction 
band and the top of the valence band contact each other 
at two points (We call them "contact points" ) in the first 
Brillouin zone [7, 8]. The Fermi energy of the electron 
system locates just at the contact point. These results 
were also supported by first-principle band calculations 

In the picture of a zero-gap system, transport phenom- 
ena of this material are understood as follows . Carrier 
density written as n oc is a characteristic feature of 2D 
zero-gap conductors with the Fermi energy located at the 
contact point [1, Q . Carrier mobility, on the other hand, 
is determined as follows. According to Mott's argument, 
the mean free path [V] of a carrier subjected to elastic 
scattering can never be shorter than the wavelength (A) 
of the carrier, so that I > A[21|. For the case that there 
scattering centers exist with high density, Z ~ A. As the 
temperature is decreased, the mean free path becomes 
long because the wavelength (A) becomes long with de- 
creasing energy of carriers. As a result, the mobility (/i) 
of carriers increases as /i cx T^^ in the 2D zero-gap sys- 
tem 1^, Q . Consequently, the resistivity per layer (sheet 
resistance Rs) is given as Rs = h/e^, which is indepen- 
dent of temperature @ . 

One problem that remains to be explained is the trans- 
port phenomena under a magnetic field [1, [l^, Ull • As 
transport phenomena of this system are extremely sensi- 
tive to the magnetic field, the samples are expected to be 
clean. Both quantum oscillations and the quantum Hall 
effect, however, were not observed in the magnetic field 
of up to 15 T. This is explained by the fact that the elec- 
tron system is not degenerated. The Fermi energy of this 
material locates at the contact point both in the absence 
and in the presence of a magnetic field. Therefore, elec- 
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trons do not have the Fermi surface. This is one of the 
characteristic features of intrinsic zero-gap conductors. 

Another important difference of zero-gap conductors 
from conventional conductors is the appearance of a Lan- 
dau level at the contact points when magnetic fields are 
applied normal to the conductive layer. This special Lan- 
dau level is called the zero-mode level. In this work, we 
succeeded in obtaining evidence of the zero-mode Lan- 
dau level in our material. The idea of this work is as 
follows: Since the energy of this level is Ep irrespective 
of field strength, the Fermi distribution function is al- 
ways 1/2. This means that half of Landau states in the 
zero mode are occupied. Note that in each Landau level, 
there are states with density proportional to B. Thus, 
the magnetic field creates mobile carriers. 

This effect can be detected in the interlayer resistance 
under a transverse magnetic field, as pointed out by Os- 
ada [13]. In this field configuration, the interaction be- 
tween the electrical current and the magnetic field is 
weak, because they are parallel. On the other hand, the 
magnetic field in this direction (interlayer) yields Landau 
carriers. Hence, the conductivity is primarily propor- 
tional to the carrier density. The effect of the magnetic 
field appears only through the change in the carrier den- 
sity. 

To demonstrate the appearance of zero-mode Landau 
carriers, we examined the inter-layer resistivity of this 
zero-gap conductor in the magnetic field. The resistance 
was measured at temperatures down to 0.06 K and at 
magnetic fields up to 7 T. 

A sample to which four electrical leads were attached 
was encased in a Teflon capsule filled with pressure 
medium (Idemitsu DN-oil 7373). The capsule was set 
in a clamp- type pressure cell made of MP35N hard alloy, 
and hydrostatic pressure of up to 1.7 GPa was applied. 
Resistance measurements were carried out using a con- 
ventional dc method with the electrical current along the 
c-crystal axis, which is normal to the 2D plane. 

Figure 1 shows the magnetic field and the angle de- 
pendences of the inter-layer magnetoresistance. Resis- 
tivity behavior is dependent on the direction of the mag- 
netic field. In the magnetic field parallel to the 2D plane 
{6 — 0°), a gradual increase in resistivity was observed. 
In the field normal to the 2D plane {9 = 90°), resis- 
tivity also increases at low fields below 0.2 T. Then, it 
decreases. Note that the resistivity drop is large: At 3 
T, the resistivity is about 1/3 of that at 0.2 T. 

To understand the data for B ± 2D plane in Fig. 1, we 
should examine zero-mode carriers quantitatively. The 
energy of Landau levels in zero-gap systems is expressed 
as -EnLL ~ ±^y2ehvp\n\B, where vp is the Fermi velocity 
and n is the Landau index (23l |. In drawing the inset in 
Fig. 1 (a), we used vp = 10^ cm/s estimated from the 
data of the temperature dependence of carrier density. 
The carrier density depends on temperature as n = aT^, 
where a = (n'^k'^/Qch'^Vp) and c is the lattice constant 
along the direction normal to the 2D plane Accord- 
ing to band calculations, however, the Dirac cone in this 
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FIG. 1: (color online), (a) Magnetic field dependence of in- 
terlayer resistance under the pressure about 1.7 GPa at 4 K. 
When B > 0.2 T, remarkable negative magnetoresistance is 
observed. As for negative magnetoresistance, fitting is done 
with Eq. (1), Rzz oc 1/(|-B| -I- Bo) (red line). Inset shows 
the Landau levels of this system, (b) Angle dependence of 
magnetoresistance. When 9 = and 9 = 180°, the direction 
of the magnetic field is parallel to the 2D plane. 9 = 90° 
is the direction normal to the 2D plane. Fitting of the data 
measured at 1, 2, and 3 T was done using Eq. (1). 



system has highly anisotropic Fermi velocity vp {4>) 0, Hi ■ 
Therefore, vp = 10''cm/s is an average of vp{4>). The 
inset in Fig. 1 (a) agrees with the realistic theory by 
Goerbig et al. [24| . 

Since the experiments were performed at temperatures 
below 10 K, the energy gap between the zero-mode level 
and the n=l level could exceed this temperature at rather 
low magnetic fields. At 1 T, for example, the energy 
gap is about 20 K. Therefore, if the experiment were 
performed at temperatures below 20 K and at magnetic 
fields above 1 T, most of the mobile carriers would be in 
the zero-mode Landau level. Such a situation is called the 
"quantum limit " . The quantum limit region is expressed 
as kBT < iJiLL- 

Carrier density in the quantum limit is expressed as 
D{B) = B/24>o, where 0o = h/e is the quantum flux. 
Factor 1/2 is the Fermi distribution function at Ep. In 
moderately strong magnetic fields, the density of carriers 
induced by the magnetic field can be very high. At 3 
T, for example, the density of zero-mode carriers will be 
10^^ m~^. This value is by about 2 orders of magnitude 
higher than the density of thermally excited carriers at 
4 K, and, in the absence of the magnetic field, it is 10^'^ 
m^^. Therefore, carrier density in the magnetic field is 
expressed as B/2(j)o except for very low fields. This large 
change in the density of zero-mode carriers is the origin 
of the negative magnetoresistance, as shown in Fig. 1. 

Recently, Osada gave an analytical formula for inter- 
layer magnetoresistance in a multilayer Dirac fermion 
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system as follows: 



A 



l^zlexp 



h\B, 



(1) 



Bn 



where A = ttH^ /2Ct'^ce^ is a parameter that is considered 
to be independent of the magnetic field if the system is 
clean. Bq is a fitting parameter, C is defined by C = 
/ po{E){—df /dE)dE using the spectral density of the 
zero-mode Landau level and po{E) satisfies / po{E)dE = 
1 [l^. Note that only lattice constant c is a material 
parameter. 

Except for narrow regions around 9 — 0° and 180°, 
this formula can be simplified to pzz = ^/{\Bz \ + Bq). 
Using this formula and assuming Bq—O.I T, we tried to 
fit the curves in Fig. 1. This simple formula reproduces 
well both the magnetic field dependence and the angle 
dependence of the magnetoresistance at magnetic fields 
above 0.5 T as shown by solid lines in Figs. 1 (a) and 
(b), which evidences the existence of zero-mode Landau 
carriers in a-(BEDT-TTF)2l3 at high pressures. 

Here, we briefly mention the origin of positive mag- 
netoresistance around 9 = 0° or 180°. In the magnetic 
field in these directions, the Lorentz force works to bend 
the carrier trajectory to the direction parallel to the 2D 
plane. It reduces the tunneling of carriers between neigh- 
boring layers so that the positive magnetoresistance is 
observed. Note that the formula (1) for 6* = 0° or 180° 
dose not correctly evaluate the effect of the Lorentz force 
and, thus, loses its validity. The value of the resistance 
peak depends weakly on the azimuthal angle. At 3 T, for 
example, the ratio of maximum value to the minimum 
value is less than 1.3. According to the calculation of in- 
terlayer magnetoresistance by Morinari, Mimura and To- 
hyama, the effect of Dirac cones with highly anisotropic 
Fermi velocity is averaged and gives rise to this small 
difference f25| . 

An apparent discrepancy of the data from the formula 
(1) is also seen at both low and high magnetic fields nor- 
mal to the 2D plane, because the model is oversimplified. 
Equation (1) was derived based on the quantum limit 
picture in which only the zero-mode Landau level is con- 
sidered. In fact, each Landau level has a finite width 
due to scattering. At a sufficiently low magnetic field, 
the zero-mode Landau level overlaps with other Landau 
levels. In such a region, the formula (1) loses its valid- 
ity. We can recognize this region in Fig. 1 below 0.2 T, 
where positive magnetoresistance is observed. This crit- 
ical magnetic field shifts to a lower field with decreasing 
temperature, as shown in Fig. 2 (a). 

The deviation of data in the high field region is much 
more serious. In this region, the resistance increases ex- 
ponentially with increasing field. This phenomenon is 
understood as follows. 

In the above discussion, we did not consider the Zee- 
man effect. The Zeeman effect, however, should be taken 
into consideration because it has a significant influence 
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FIG. 2: (a) Magnetic field dependence of interlayer resistance 
for the magnetic field normal to the 2D plane at several tem- 
peratures from 0.06 K to 9.6 K. (b) Temperature dependence 
of B\ when resistance is assumed to be Rzz oc exp(_B/_Bi). 
(c) Temperature dependence of the width of Landau levels 
(2fi/r). 



on the transport phenomena at low temperatures. In 
the presence of a magnetic field, each Landau level is 
split into two levels with energies EnLL ± Ai?, where 
AE—pbB is the Zeeman energy. This change in the en- 
ergy structure gives rise to a change in the carrier density 
in Landau levels. In particular, the influence on the zero- 
mode carrier density is the strongest, because the energy 
level is shifted from the position of the Fermi energy. 
The value of the Fermi distribution function varies from 
fiEp) = 1/2 to fiAE) = l/ieM^E/kBT) + 1). At 
low temperatures where /cb^ < AE, this effect becomes 
important. It works to reduce the density of zero-mode 
carriers and, thus, increases the resistance. 

The Zeeman energy when _B = 1 T is about 1 K. There- 
fore, in the experiment performed at 1 K, the deviation 
of experimental results from Eq. (1) is expected to start 
around 1 T. This is confirmed in Fig. 2 (a). At 1.8 K, 
for example, the deviation is prominent in fields above 2 
T. This critical field shifts to about 0.3 T at 0.06 K. 

A definite evidence showing that the anomalous in- 
crease in the resistance at high field in Fig. 2 (a) is due to 
the Zeeman effect is given by examining the slopes of the 
curves at high fields. In this region, the magnetic field de- 
pendence of resistance is expressed as Rzz oc exp{B/ Bi), 
where Bi is a parameter that depends on temperature. 
Bi for T =1.8 K is estimated to be about 2.7 T. At 0.06 
K, it decreases down to about 1.1 T. In Fig. 2 (b), we 
plot the temperature dependence of Bi. Above 1 K, the 
T vs Bi curve is close to a line Bi — k^T/pB- This is 
strong evidence that the Zeeman splitting of zero-mode 
Landau levels is the origin of the resistance that obeys the 
exponential law as Rzz oc exp(2/iB^/2A:BT). In contrast, 
the behavior of Bi below 1 K is understood in terms of 
the width of Landau levels. According to the theory of 
Osada, p-rB / (h/r) ~ 1 at the magnetoresistance mini- 
mum [22| . By using this relationship, the energy width 
of zero-mode Landau levels (2/i/r) is estimated as shown 
in the inset in Fig. 2 (c). With decreasing temperature. 
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the energy width decreases hnearly down to about 0.7 
K, where it almost saturated. At 0.06 K, it is estimated 
to be about 0.5 K. This energy is much higher than the 
thermal energy. This is the Dingle temperature in the 
present system. 

In conclusion, we provided ample evidence that a- 
(BEDT-TTF)2l3 under high hydrostatic pressure is an 
intrinsic zero-gap conductor with the Dirac-type energy 
band. This is the first work that investigated the effect of 
the zero-mode Landau level on the transport phenomena 
of bulk zero-gap conductors. We succeeded in detecting 
the zero-mode Landau level that appeared at the contact 
points of Dirac cones under the magnetic field normal 
to the 2D plane. The large negative magnetoresistance 
observed in the interlayer resistance in longitudinal mag- 
netic fields is due to the increase of the zero-mode car- 
riers. Good agreement between experiment and theory 



was obtained. At high magnetic field, spin splitting of the 
zero-mode level becomes important, and the resistance 
obeys the exponential law. This is the characteristic fea- 
ture of transport in multilayer zero-gap conductors. In 
graphite systems, a monolayer sample (graphene) is in- 
evitable to realize a zero-gap state. On the other hand, in 
the present work, we demonstrate that bulk crystals can 
be 2D zero-gap systems. a-(BEDT-TTF)2l3 provides us 
a suitable testing ground for the transport of the multi- 
layer massless Dirac fermions system. 
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